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ABSTRACT 


Exhaustive studies have been conducted in the field of 
load flow analysis. The results of these studies have pointed 
out the advantages of the fast converging methods like 
Newton-Raphson especially in polar coordinates and Fast Decoup- 
led method . The former method has the disadvantage of large 
memory requirement and greater computation time. These 
methods (in particular Newton-Raphson) become a practical 
tool only when sparsity of the coefficient matrix and bus 
admittance matrix are exploited. The use of ordered 
elimination reduces the computation time further. 

In this thesis, programs have been developed to 
perform LF studies by all the three methods such as, , 

D.C. and FD. Sparsity ordered eliminations are the key feature 
of the program. A comparative study of these three methods 
with reference to computation time and memory requirements 
is also given. 



CHAPTER 1 
INTRODUCTION 


The modern trend is to form a grid system of all the 
available energy sources i.e. towards interconnecting all types 
of generating stations. This provides the greatest advantage 
of meeting the load supply demands economically at all times. 

The power supply undertakings must keep pace vjith the load 
growth. In addition to the above, care must be taken so as 
not to overload the interconnecting systems resulting into their 
instability. 

Load flovj analysis is very important Vi/hen new components 
or additions to existing ones are considered. With proper 
and accurate load flow studies, the interruption of power can be 
minimised. Load flow calculations are necessary at the initial 
stage for the purpose of planning, operation and control. It 
provides voltage magnitude and phase angle at each bus and 
power flows including line losses in each element of the power 
system network. Apart from determining the steady state 
operating conditions of a power system network for the purpose 
of planning, operation and control, load flow calculations also 
provide initial conditions for transient stability studies. 

Prior to the advent of digital computer, load flow 
studies were perfoimed on A.C. calculating boards i.e. network 
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analysers, A calculating board is a single phase scaled 
down model of a balanced three phase system. The board being 
made up of a number of elements viz. resistances, inductances 
and capacitances, all of which are adjustable, along with a 
number of sources and measuring instruments. Initial adjust- 
ments in this case usually take a lot of time since each 
adjustment at any bus affects values of pertinent quantities 
at other buses. In addition to this, considerable amount of 
time is lost in recording observations.. 

The appearance of the digital computers revolutionized 
the whole concept of load flow calculations. Mathematical 
model (i.e. equations) which were once thought to be cumbersome 
and of purely theoretical interest became practically feasible. 
The ease with which computers can handle arithmetic 
operations gave a boost to the numerical methods. The mathe- 
matical model for the purpose of load flow studies is a set of 
non-linear algebraic equations. The non linearity of the 
system of equations defies an exact analytical solution and one 
must resort to some iterative techniques which will render a 
sufficiently accurate numerical solution. There is no dearth 
of numerical techniques available, only the enormous computa- 
tional effort is a deterrent, but with the coming of the 
digital computer, it is no longer a stumbling block, for now 
the problem is to develop an algorithm for solving these 
equations on the computer. 
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The first practical methods to solve these power system 
network equations on a digital computer, appeared in literature 
in 1956 [ 1 , 2 ]. These methods (one of the methods was the 
gauss-seidel technique) required minimum storage and hence 
were well suited to the first generation computers. However 
these methods were slow in convergence and thus not very well 
suited to handle large systems. Any method v;hich has to 
handle a large system must possess the following two key 
features, 

1. Nominal storage requirements 

2. Reliable and fast in convergence. 

The Newton-Raphson method's quadratic convergence 
property was highlighted around the same time [3,4] but was 
found to be computationally uncompetitive. The application of 
sparsity programmed ordered elimination by Tinmy and Walker 
to the Newton-Raphson method reduced the storage requirement and 
also optimized the computation time to such an extent that 
Newton-Raphson method gained popularity over and above other 
methods [5], and has now come to be widely regarded as the 
general purpose load flow approach [6j. The decoupled and 
fast decoupled load flow techniques are modifications of the 
Newton-Raphson method which exploit the loose physical 
interaction between MW and MVAR. flows in a power system. 

Storage and computation time are further minimized in the 
above mentioned methods, without appreciable loss in accuracy. 
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Present day power systems are large and complex because 
of greater interconnection. To analyse such a large scale 
system on a digital computer with limited memory application 
of sparsity oriented ordered elimination techniques are needed* 

Keeping the above factors in view, programs are 
developed for the three methods viz. Nowton-Raphson, Decoupled 
and Fast Decoupled. These programs have been tested for a 
100 bus 128 line system of UPSEB. Programs are capable of 
handling a larger system; data storage requirements of the 
large system are the limiting factors which dictate system size 
that can be simulated on a particular digital computer. 

Although the use of magnetic tapes can overcome this problem 
to some extent, one has to pay in terms of speed. The main 
features of the programs developed in this thesis are: 

1, User oriented input/output format 

2. Storage of only non-zero elements of 

3, Storage of only non-zero elements of Jacobian 

4. Ordered elimination of the Jacobian equation 

The chapter-wise summary of the work covered in this 
thesis is given as follows. 

Chapter 2 is devoted to the theoretical aspects of the 
methods used viz. Newton-Raphson, Decoupled and the Fast 
Decoupled method. A brief account of each method and their 
relative merits are also discussed in this chapter. 
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Chapter 3 deals with sparsity . ordered 
elimination. A general description of the technique and in 
specific, its application to power system problem has been 
given. 

Chapter 4 deals with the case study of the following 

systems 

14 bus 20 lines IEEE system 

57 bus 80 lines IEEE system 
and 100 bus 128 lines UPSEB system. 

The advantage of sparsity ordered elimination are 
elaborated by compare sion of results for the three systems in 
relation to memory requirements and computer time. Detailed 
flow chart for all the methods used as well as results etc. 
are given. 

Chapter 5 concludes with the specific findings in this 
thesis along with future scope of the work. 



CHAPTER 2 


LOAD FLOW TECHNIQUES 

2.1 INTRODUCTION: 

This chapter deals with the currently favoured methods 
for load flow studies. A literature survey will reveal a host 
of algorithms which have been suggested from time to time to 
solve this problem of load flow analysis. An excellent review 
of the major portion of work done in this field has been 
given in [7j, In general it is difficult to point out the 
best method fof a particular application. The relative 
properties and performances of different load flow methods can 
be influenced substantially by the types and size of the 
problems to be handled and also by the computing facilities 
available. Any final choile is invariably a compromise between 
the various criteria of goodness by which the load flow methods 
are to be compared with each other. Every such criteria is 
directly or indirectly associated with financial cost. This 
chapter spells out the details of the load flow problem and 
the numerical techniques for its solution, 

2.2 LOAD FLOW STUDIES; 

The objective of the load flow study is to 

determine the phase angle and reactive power on each P-V bus 
and the phase angle and voltage magnitude at each P-Q bus 
subject to the constraints on the real and reactive power at 
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P-Q buses and the real power and voltage magnitude at the 
P-V buses. Based upon this it is possible to classify the 
buses into three categories, 

2.3 BUS CATEGORIZATION: 

The buses are categorized depending on the quantity 
specified at the bus 

a) Load or a P-Q bus 

b) Voltage controlled or a P-V bus 

c) Slack or swing bus 

a) Load or a P-Q bus : For this type of a bus, we know a priori 

Pj and Qj and specify P^ and . In effect we thus 

specify the bus injections P^ and Q^. Solution of the load flow 

equations will render jv^^j and A load bus which due to its 

lack of generating equipment, is characterized by zero P- and 

^i 

Qq^ evidently falls in this category, 

b) P-V or a voltage controlled bus : For this type of a bus we 

know a priori P, and Q, and specify jV-j and P^ . In effect, 
we thus specify the bus powers P^^. Solution of the load flow 
equations render Q^ (and hence Qq_ ) and This is called 

a voltage controlled bus because its voltage can be controlled, 

c) Slack or swing bus : This is the reference bus where the 
voltage magnitude and phase angle are specified. One of the 
generator with the maximum real power capabilities must be 
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selected as the swing bus to provide for the additional real 
and reactive power to supply line losses because these are 
unknown till the final load flow solution is obtained. The 
variables of interest at this bus are the real and reactive 
power. 

Assuming balanced 3-phase conditions, which is usually 
done for the purpose of load flow studies, the transmission 
system can be represented by its positive sequence network. 
The nodal admittance matrix can be expressed as follows 


^BUS “ ^BUS ^BUS 

(nxl) (nxn) (nxl) 


( 2 . 1 ) 


The above equation can be w/ritten in the following 
form for a node, 

n 


I 


P 


Z 

q=l 


Y V 

pq q 


( 2 . 2 ) 


p = 1,2, ...,n 


This equation simply states that the currents at any 
node or bus is the algebraic sum of all the currents entering 
or leaving the node. The power at any bus is calculated by 

the VI* product. 

P P 




n 

E 

q=l 




(2.3) 
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separating equation (2.3) into the real and imaginary parts 
gives us the expressions for real and reactive powers i.e. 


(2.4) 


Q 


■pc( 


[v^ 

n 

£ 


v; 

] 

p 

q=l 

n 

pq 

q 


E 


v* 

] 

p 

q=l 

pq 

q 


following 

'pq ^'^pq 

substitutions for Y _ 

pq 


(2>5) 


V. 


jVpt (cos 0p + j sine 9^ ) 


Vq = jV^I (cos 9^ + j sine ) 


equations (2.4) and (2.5) become 


n 


V„1 2 ((G^ cos 9 + B sine 9 )lv 1) (2.6) 


q=l 


pq 


pq pq 


pq ' q' 


n 


IVnl ((Gpq =1"® «= e„„)|v„i) (2.7) 


q=l 


pq pq 


pq'* q' 


Let US examine the number of knowns and unknowns at 
the three type of buses. 


suppose total number of buses 
slack bus 

Number of P-Q buses 
Number of P-V buses 


= N 
= 1 
= M 
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If V and © are knomat all the buses we can find out P 
and Q at all buses using equations (2.6) and (2.7) i.e. V and 
9 are the state variables. 

Total number of possible unknowns = 2N 
As voltages at all P-V buses are known and also 
at slack bus V is assumed 1,0 p.u. and angle 0^ = 0*^, 

Number of knowns = N-M--1+2 

= N~M+1 

Number of unknowns = 2N-(N-M4-l) 

= N+M-1 

Hence (N+M~l) equations are needed to solve for the 
unknowns. For each load bus P and Q are known so we can write 
two equations at each P-Q bus. Also P is known at each P-V 
bus so one equation for P can be written for each P-V bus. 

Number of equations for (N-M-1) P-V buse« « N-M-1 

Number of equations for M P-Q buses = 2M 

Total number of equations = 2M+N-M-1 

= N+M-1 

Thus the number of equations is equal to the number of 
unknowns [Note this has been possible, only if Q i.e. 

(0 - © ) is treated as one unknown by taking one of the 

P Q 

buses as reference]. 
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Bus constraint equations are 


^Pp 



pcal 

P 


(2.9) 


AQp = QpP-Qp®^ (2.10) 

where superscript ' sp* and ’cal’ stand for specified and 

^ ^ If 

calculated respectively. and are obtained from 

the equations (2,6) and (2.7). As can be seen by the 

appearence of cos 9 and sine 0 terms in the expressions 

pq pq 

C 3 JL C 3 ^ 

for Pp and , it is a system of non-linear equations 

and one has to resort to numerical techniques to obtain a 
solution. The solution of these equations for V’s and 0’s 
is the load flow problem. 


2.4 NEWTON RAPHSON METHOD: 


V'llhen there is no mismatch between the specified and 
calculated powers equations (2.9) and (2.10) [in matrix 
notation] become 


A P 
AQ 


= 0 


( 2 . 11 ) 


Applying Newton-Raphson method we have 





( 2 . 12 ) 
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In short the above can be written in the form 


H N 


u T Av I 

M L jq-; 


(2.13) 


where 


H N 


M L 


is called the Jacobian matrix 


H — Partial derivatives of P w.r.t. 0's 

L — Partial derivatives of Q w.r.t. Vs 

N — Partial derivatives of P w.r.t, V s 

M — Partial derivatives of Q w.r.t. 0's 

The A V‘s are divided by jvj and corresponding elements 
of Jacobian are multiplied by 'V to bring about a symmetry in 
the elements of the Jacobian. 

It can be shown that 
for p j4q 

Hpq = Lpq = I’^ql <°pq ®pq ’ ®pq ®pq^ 

% = -«pq= IVpl IVql Spq + ®pq ®Pq> 


For p = q we have 

Hpp = -Qp - Bpp I'^pl 


(2.16) 


Q - B V 
> pp ' p' 


(2.17) 
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pp 

= P^ + 

G 

V ■ 

p 

pp 

' p* 

‘pp 

= p ~ 

G 

Iv 1' 

p 

pp 

p 


Where 

P and Q 

H j 


(2.6) and (2.7). 


calculated from equations 


(2.18) 

(2.19) 


The solution of equation (2.13) gives us the A 0‘s 
and A V’s which are used to update earlier estimates of 0’s 
and V’s and the process is repeated till the mismatch A P 
and AQ become less than a pre-assigned tolerance value 
When this is achieved, the iterative process is stopped as the 
desired accuracy has been obtained. 

2.5 STEPS FOR THE NEWTON-RAPHSON ITERATIVE SCHEME; 

1. If nothing is available regarding the actual value of 
variables at the buses assume a flat start, assign V at 
all buses equal to slack bus voltage and angles equal 
to slack bus angle i.e. zero. Set iteration count *K’ 
to one. 

2. Calculate and (using equation (2.6) and (2.7)) 

with values of V’s and @’s as in step (l). 

3. Calculate power mismatch at all buses using equation 
(2.9) and (2.10). 

4. Test for convergence by checking power mismatch. If 

A P’s and AQ’s at all buses are less than a pre-defined 
value ^ ,we jump out of the iterative loop and go to 
Step ( 10 ) . 



15 


5, Check if the number of iterations has exceeded the 
predefined value ’UMAX’ (say), if it has exceeded go 
to step (12). 


6. Calculate the elements of the Jacobian using equations 


(2*14)through (2.19). 

7. Solve equation (2.13) for 0‘s and 
the direct methods of solution (e.g. 


, - X 

* s using one of 
Gaussian elimination) 


8. Update the voltages and angles at all the buses using the 
correction factors obtained in step (7). Increment 
iteration count by ‘1’ 

jvjK+l ^ j^jK ^ Ivj (2,20) 

0^*^^ = 0^ + A 0^ (2.21) 

9. With the voltages and angles as given equations (2,20) 

■end (2.21) start the (K+1) iteration i.e, go to 

step (2). 

10. Using the latest voltage and estimates , calculate slack 
bus power, line flows and line losses. 

11. Go to Step 13 . 

12. Convergence not obtained in ’K’ iterations. 

13. Convergence obtained in 'K' iterations. Print bus status^ 
line flows, line losses. 
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The main disadvantage of this method is that the storage 
requirements and computation work involved is enormous. For a 
'N’ bus system with 'M* P-Q buses the order of Jacobian is 
(N+M-l). Thus for a typical 100 bus problem with 19 P-V buses 
including slack bus, which has been carried out in this 
thesis, we require 32.4 K of computer memory for storing the 
Jacobian matrix. Storage of data, bus admittance matrix etc. 
are over and above this. Bus. admittance matrix for a 100 bus 
system will contribute towards a storage requirement of lOK. 

Thus Jacobian and bus admittance matrix together take the 
major portion of total storage requirements for any problem. 

With full storage schemes the solution is limited to small 
problems because of memory restrictions. The Newton-R-aphson 
method together with sparsity and ordered elimination technique 
[5] is a powerful tool for obtaining load flow solution, as 
it optimizes memory requirement as well as computation time. The 
number of iterations required for solution is virtually 
independent of problem size. This is strictly true for 

programs without additional features like automatic tap adjust- 
ment of a transformer, Q limit checks etc. which may require 

additional iterations. A program adjusted for Q limits may 
take an additional two or more iterations. 

2.6 DECXiUPLED METHOD; 

In all the decoupled methods the load flow equations 
have been derived from the Newton-Haphson formulation in pola; 
coordinates to reduce memory requirement and computational 
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efforts. These methods are based on neglecting the coupling 
terms M and N of the Jacobian matrix in the Newton-Raphson 
method, on the assumption that the coupling between real bus 
power versus bus voltage magnitude and reactive power versus 
bus voltage angle is relatively weak. Any such approximations 
to the Jacobian inevitably sacrifices the true quadratic 
convergence property, but compensating computational benefits 
can accrue. Based upon these assumptions equation (2.13) 
reduces to two sets of independent equations for P’s and Q’ s. 

CAP ] = H [/\Q] (2.22) 

[AQ ] = L [^] (2,23) 

Equations (2,22) and (2,23) are formulated and solved 

successively. The latest values of @ are used to solve for V. 

The decoupled method converges as reliably as the formal 

Newton-Raphson Method, although it takes more number of 

iterations to achieve accuracies compatable to the Newton’s 

method. This however is not necessary as convergence to 

practical accuracies takes more or less the same number of 

of 

iterations. The saving in terms £ memory requirements is 
nearly 75/< for Jacobian element storage although overall 
saving of the memory is only of the order of 40-50;/i. The 
computation time per iteration is also 10-20°A less than 
Newton-Raphson Method. 
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2.7 FAST DECOUPLED METHOD: 


The decoupled method can be further simplified without 
appreciable loss of accuracy [7, 8]. In practical power system 
the following assumptions hold good. 


1 . 

2 . 

3. 


0 


pq 


is small . 


G sine © << B 

pq pq pq* 


«p « Bpp |V 


Applying these assumptions to equations (2,22) and 
(2.23) [reproduced below]. 




[Ap] = 

H 




[Aq] = 



We have 





[A 

P] 

= [v B'V] 

][Ae] 

(2.24) 

[A 

Q] 

- [V B”V] 


(2.25) 


The elements of the matrix B' and B‘ ' are strictly 
elements of [-Bj, The decoupling process is given a final 
shape by. 

(a) Omitting from [B’ ] the representation of those network 
elements that predominantly affect MVAR flows i.e, 
shunt reactances and off nominal in phase taps. 



(b) Omitting from [B''] the angle shifting effects of phase 
elements# 

(c) While calculating for bus taking the left hand ‘V’ 
terms (for bus) in equations (2.24) and (2,25) on to 
the left hand side of the equations and then in equation 
(2.24) removing the influence of MVAR flows on the calcu 
lations of A © by setting all right hand 'V* terms to 

1 p.u. 

With these assumptions the relevant equations for 
Fast-Decoupled load flow are 

= [B'] [Ae] (2.26) 

= [B"] [Avj (2.27) 

This method though not possessing the true quadratic 
convergence of the Newton-Raphson method^ converges very fast 
as the time per iteration is very less. It is as reliable as 
the Newton-Raphson method within the acceptable limits of 
accuracy. Adjusted solutions, to incorporate all other 

additional features, in this case, will take more number of 

Z ince 

time per iteration is very less compared 
to Newton-Fiaphson method the overall computation time is 
not affected significantly. 
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In this chapter we have outlined the various methods 
of current interest. The methods in themselves are not new 
but form a powerful tool when sparsity of the Jacobian 
matrix is exploited. Implemented as such^ they may not be 
able to handle systems of 500 bus or more (especially 
Newton-Raphson method) whereas using sparsity, we can handle 
system sizes of 1000 buses and above with little difficulty. 



CHAPTER 3 


SPARSITY AND OPTIMAL ORDERIN3 
3.1 INTRODUCTION: 

The sparsity occurs in some form in most of the 
physical systems such as communication network, current 
theory , family trees, organization structure and sociograms. 
Let a physical system be described by a set of 'n' algebraic 
linear equations of the form 

[a] X = y (3.1) 

The problem is to determine the solution vector x by 
Gaussian elimination method such that the . computational 
efforts and hence, the time of computation i.e. the cost is 
minimized. Following Von Neuman, the number of multiplication 
required to obtain solution is counted as a measure of 
computing time. Therefore if only the number of multiplica tio 
is to be counted, a reduced matrix ‘M* of the coefficient 
matrix A [Eqn. 3.1 ], whose elements are defined as 


=1 if m. . ^ 0 


(3.2) 

contains all the required information for solving the problem. 


0 if m. . = 0 


"i” W 

Let the number of multiplication to process the i 
row be m^^ and therefore, for the entire system, the total 
number of multiplication 
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n 

0 = S m. (3.3) 

i=l ^ 

where n is the number of equations i.e. order of the system. 


3.2 FUNCTIONAL EQUATIONS: 


Optimal elimination is actually a topological problem 
which can be formulated using notation from graph theory. Some 
systems such as electrical networks may be thought of being 
their own graph, thus a picture of one of these systems with 
slight modification could serve as its own graph inspite of 
the fact that there may be more than one scaler quantity 
associated with each node, other systems such as those arise 
from difference equations may have no direct graph. For these 
systems, the following procedure is adopted to construct its 
graph. With each equation in the coefficient matrix *A’ 

(eqn. 3.1), there is associated a node in the system graph 
and with each non-zero term, 


ay for 1 = 1, ... n 
j = ;., . . . n 


(3.4) 


there is associated an undirected branch between the ith and 
the jth node. 

The system graph will be referred to as ’G*j during 
the elimination process it is modified, just as the rows of 
the coefficient matrix 'A’ arc modified. Let is a graph 
obtained by eliminating the ith node from system graph G 
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[i.e, processing ith row of the corresponding coefficient 
matrix 'A’J. Let V1(G) be defined to be the number of 
multiplication required to solve optimally the system [i.e, 

[a]x = y] whose graph is ' G’ and let ’(i) be one plus the degree 
of i"^^ node in the graph ' ^G’ . From this, it is clear, that, 
W(G) is a minimum value of 0 and e^ is the number of 
multiplications required to eliminate itg row from the given 
system whose graph is G, Then, 

W(G) = Min 0 (3.5) 

W(G) = Miin [e^ + VJ(G^) ] (3.6) 

where is the permutation of ordering. 

Bellman uses the term 'policy* to describe a specific 
permutation i.e. a certain policy results in a permutation for 
which it is then possible to evaluate the number of multipli- 
cations or work. The optimal policy corresponds to the 
minimum work. 

Following Bellman, it is possible to choose any initial 
policy i.e. method of ordering the nodes i.e. equations or 
rows of corresponding coefficient matrix and proceed iteratively 
to obtain the solution of the above equations i.e, eqn. (3.6) 
whose solution is unique even though the optimal policy may 
not bo unique. 
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Let W^[G] be the number of multiplications needed 
using initial policy* Then we have from the equation (3.6) 

Wj^[G] = Min [e^ + Wj^_^[G^]j (3.7) 

i 

Here Wj^Cg] is the number of multiplications required to 
solve optimally the system having 'n' equations i.e, whose 
coefficient matrix has ‘n* rows and G is the corresponding 
graph. The solution of these equation which is dynamic 
programming will yield the following result. At each step in 
the elimination scheme, eliminate that node next which has 
the smallest degree. 

Such problems can easily be formulated and solved by 
the principle of dynamic programming which is developed by 
Richard Bellman because these problems belong to a category 
known as the multistage decision process, typical example 
being that^travelling sales man problem. Here we take the 
initial decision which is arbitrary and based upon this 
decision all other decisions are optimal, say in this 
particular case, the initial decision is that the ith row of 
the coefficient matrix ’A’ is processed first i.e,, ith 
node of the corresponding graph 'G* is eliminated first; 
because of taking this decision, the cost involved e^^ where 
e^ indicates the number of multiplications needed to process 
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4* H 

the row of coefficient matrix A, will be the measure of 

the cost to process the ith row. Because of taking this 
decision, the graph ' G’ will change to and number of nodes 
will become (N-l) and hence the formulation using dynamic 
programming will yield the result, 

Wj^(G) = Min (ej_ + Wj^_^(G^)) (3.8) 

i 

N = 1 ... n 

The main advantage of using this formulation is; at 
any stage we deal with only one variable i.e. instead of 
solving all the n variables together, they are solved one at 
a time, however n number of times. 

3.3 DIRECT SOLUTION OF SPARSE NETWORK EQUATIONS BY OPTIMALLY 
ORDERED ELIMINATION; 

For the sparse systems, which normally occur in power 
system network formulation, solution is obtained by optimally 
ordered elimination. This method consists of two parts 
[ 9 , 10 , 11 ], ■ 

1) A scheme of recording the operation of triangular decompo- 
sition of a matrix such that repeated direct solution can 

be obtained without repeating the triangularization process. 

2) A scheme of ordering the operation such that it tends to 
conserve sparsity of the original system. 
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The first part of the method is applicable to any 
matrix. However the application of the second part i.e, 
ordering to conserve sparsity is limited to sparse matrix in 
which the pattern of non-zero elements is symmetric and for 
which an arbitrary order of decomposition does not affect 
adversely the numerical accuracy, such matrices are normally 
characterized by a strong diagonal. The coefficient matrix 
in the case of the load flow problem belong to this category 
where more than 90'A elements are zero at off diagonal 
locations. Let us take the equation 


[a] X = y which can bo expanded as 
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The matrix 'A* is changed to augumented matrix A by 
adding in (n4-l)th column, the known constants of column 
vector Y. By factored LU decomposition of the coefficient 
matrix, we obtain the following matrix known as the table 


of factors 
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where the elements of the matrix are defined belov; 
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ij 
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ij 


Wheri the matrix to be decomposed is sparse the order in 
which the rows are processed affect the number of non-zero 
terms in the upper triangular matrix. If a programming scheme 
is such that it processes and stores, only the non-zero terms, 
a great Swing in operation and memory can be achieved by 
keeping the table of factors as sparse as possible. The 
absolute optimal ordering scheme would result in the least 
terms in the table of factors. 


However the absolute scheme of ordering has not been 
developed as yet, we give below the following effective 
scheme of near optimal ordering. 
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1. In this scheme the coefficient matrix of a physical 
system is ordered before hand. Here the rows with only one 
non-zero element at the off diagonal locations is numbered 
first-row with two non-zero elements is numbered two 

and so on. Finally the row with the maximum non— zero elements 
is numbered last. The rows of coefficient matrix A in the 
process of elimination, are processed in this sequence. From 
the graph point of view, a node with a degree one is numbered 
one, a node with a degree two is numbered two and finally the 
row with the highest degree is numbered last. This algorithm 
is simple to program and fast to execute , however the main 
disadvantage of the algorithm is that it does not take into 
account the changes in the pattern of non-zero elements in 
the coefficient matrix, 

2, This algorithm has been derived by using the technique of 
dynamic programming by R. Bellman, In this algorithm, in the 
process of elimination, we eliminate that row next which has 
the minimum number of non-zero elements in the off-diagonal 
locations. From the graph point of view, we eliminate that 
node next which has minimum degree. This algorithm, even 
though, being more complex than the first one, is certainly 
more efficient because it takes into account the changes in 
the pattern of non-zero elements in the process of elimination. 
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3. In this algorithm, in the process of elimination, 
eliminate that row next whose elimination will introduce 
minimum number of non-zero elements in the off diagonal 
locations. From the graph point of view in the process of 
elimination, eliminate that node next whose elimination will 
introduce minimum number of new links in the system graph. 

This algorithm has not been used by us because it takes more 
time compared to (2) . However, if the criteria is only to 
optimize the computer memory with cost having no consideration, 
this is certainly the best. 

Algorithm (2) which claims to optimize both the 
almost 

computer memory and^ the computer time has been used 

by us. The input information in this case is a list by rows 
of the column numbers counting off diagonal non-zero terms 
(i.e. branches). This scheme no doubt is more efficient 
than the first one. 



CHAPTER 4 


LOAD FLOW ANALYSIS: CASE STUDIES 

4.1 INTRODUCTION: 

This chapter presents the load flow studies for the 
following systems. 

1. 14 bus 20 lines IEEE system 

2. 57 bus 80 linffi IEEE system 

3. 100 bus 128 lin^ UPSEB systems 

These systems have been studied using the following 
methods . 

1. Newton-Raphson method in polar coordinates 

2. Decoupled method in polar coordinates 

3. Fast Decoupled method 

The choice of a particular method invariably depends 
upon the following factors, 

1. Memory requirement 

2. Speed 

3. Accuracy 

.and 4, Convergence criterion 

An attempt has been made in this chapter to compare 
the three methods based upon above mentioned criterion. The 
results of the systems studied and their significance 
are also discussed. The details of the study have been 
categorized method-wise. 
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Memory requirement and computer time invariably 
dictate the choice of method for load flow studies i.e, 
why NR method in polar coordinates has been chosen. 

Accuracy and quadratic convergence properties of this 
method are offset by the memory and computational requirement 
Although programming technique is important in all load 
flow methods for obtaining fast execution and economy in 
storage, it is the cornerstone of methods such as Newton - 
Raphson. Thus in the: case sparsity oriented programming 
makes all the difference, for without efficient storage and 
execution this method loses all its charm. To 

emphasize on the importance of sparsity oriented programming 
for these methods (especially NR method) two sets of 
programs are developed. 

SET I : Full storage mode and gaussian elimination for 
solving the load flow equations. 

SET II ; Storing only non-zero elements of the Jacobian and 
Bus admittance matrix and ordered elimination of 
the load flow equations. 

Each of the above mentioned sets offers a choice of 
three methods viz Newton-Raphson, Decoupled and Fast Decoup 
The details of memory requirement and computation time- for 
method [for the systems studied] with and without sparsity 
oriented programming are given. Each method will be taken 

CENTR ’ "ARY 

/ ; ■ ^ ■■ 
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up and studied with reference to the four factors mentioned 
before. The results of the sample systems are used as a 
means of comparing various criteria’s, 

4.2 NEVilTON-RAPHSON METHOD : 

The three systems are solved using this method. The Compi 
tation time for different systems (with and without the use 
of sparsity oriented programming) are listed in Table 4,1, 

If we consider the 14 bus system and compare the per iteration 
time in Case I and Case II we find that the difference does 
not justify the extra efforts sparsity oriented 

programming, but a glance at the results for 57 bus and 
100 bus system will speak otherwise. The iteration time in 
Case I is roughly five times that of Case II for a 57 bus system 
and twenty five times for a 100 bus system respectively. 

The saving in terms of memory requirement is also 
tremendous. Table 4,2 gives the memory saved with sparsity 
oriented approach, [Only Jacobian and bus admittance matrix 
requirements are compared as they take the bulk of storage 
space. The data storage requirements being same for both 
cases] . 

The memory requirement and computation time per iteration 
for case II increases linearly with the number of buses. 

In contrast to this, for case Id the memory requirement is 
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Table 4.1 




14 bus 

system 



Type 

C.P.U. 

Time 

No . of 
Iterations 

Time per 
Iteration 

Specified 

Tolerance 

Achieved 

Tolerance 

Without 
Sparsity 
(Case l) 

1.19 

3 

0.396 

0.001 

0.00011 

With 

Sparsity 
(Case II) 

0.95 

3 

0.316 

0.001 

0.00011 



57 bus 

system 



Type 

' C.P.U. 
Time 

' No. of 
Iterations 

Time per 
Iteration 

1 Specified 
Tolerance 

Achieved 

Tolerance 

Without 
Sparsity 
(Case I) 

34.02 

4 

8.505 

0.001 

0.00015 

With 
Sparsity 
('^ase II) 

6.98 

4 

1.745 

1 

0.001 

0.00015 


100 bus system 


Type 

C.P.U. 

Time 

No . of 
Iterations 

Time per 
Iteration 

Specified 

Tolerance 

Achieved 

Tolerance 

Without 
Sparsity 
(Case I) 

624.84 

7 

89.26 

0.001 

0.00015 

With 

Sparsity 
(Case II 

25,80 

) 

7 

3.69 

0.001 

0.00015 




















Table 4.2 



• 

c 

o 

•r-l 

-p 

O 

M 

O 

c 

•H 

c 

•H 

X 

Q> 

T5 

C 

•H 


M 

O 

4-1 

T5 

CD 

ID 

CD 

O 

c 

CD 

O 

•4-> 

05 

CD 

jC 

4^ 

CO 

CD 

X3 

D 

r-H 

O 

c 

•H 

CO 

•Ht 

x: 

H 

* 


34 


Table 4.3 


1 4 bus system 


Type 

C.P.U. 

Time 

Iterations 

Time per 
Iteration 

Specified 
^ S 

Achieved 

Without 
Sparsity 
( Ca se I ) 

3.0 

15 

0.2 

0.001 

0.00097 

With 

Sparsity 
(Case II) 

2.16 

j 

15 

0.144 

0.001 

0.00097 


57 bus system 


Type 

C.P.U. 

T ime 

Iterations 

Time per 
Iteration 

, — 

Specified 

G S 

Achieved 

6 A 

Without 
Sparsity 
(Case I) 

14.56 

7 

2.08 

0.06 

0.051 

With 

Sparsity 
(Case II) 

10.38 

7 

1.48 

0.06 

0.051 


100 bus system 


Type 

C.P.U. 

Time 

Iterations 

Time per 
Iteration 

Specified 

d S 

Achieved 

6 A 

V'iithout 
Sparsity 
(Case l) 

75.54 

9 

8.4 

0.001 

0.00048 

i 

With 

Sparsity 
(Case II) 

32.0 

9 

3.56 

0.001 

0.00048 






















Table 



3 ( 


37 


2 

roughly 5N , where N is the number of buses. From the 
Table 4.1 and 4.2, it can be inferred that the Newton-Raphson 
method realizes its full potential only when it is used with 
sparsity ordered elimination, especially for a large scale 
system. 

4.3 DECOUPLED METHOD: 

Table 4.3 gives the computation time for the three systt 
It is interesting to note that the decoupled technique, 
saves substantial amount of time, compared to Newton-Raphson 
method, for Case I, (compare tables 4.1 and 4.3 for Case I) 
but this saving is almost negligible when we compare for 
Case II, e.g. for 100 bus system, the iteration time for the 
Newton-Raphson method (case II) is 3.69 seconds while the 
corresponding time for Decoupled technique (case II) is 
3.56 seconds. Table 4.'4 gives the memory requirements for the 
Decoupled method. A comparison of Table 4,2 and Table 4.4, 
reveals the following. 

For Case I ; 100 bus system , Decoupled method requires 

29801 words of memory as compared to 52400 in 
Newton-R-aphson method. This amounts to a saving 
of 50/.. 

For case II ; 100 bus systomj Decoupled method requires 2473 

words as compared to 5126 in Newton-Raphson 
method. 


X 
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Although, the saving in this case is of the order of 
about 40 /. , its significance is not much because the 
absolute memory requirement has come down to a low level 
because of sparsity. Thus, when the sparsity and ordered 
eliminates are used, the Decoupled method is ruled out 
because memory requirement and computer time are almost 
same for both methods, hence one vi/ould prefer to make use of 
the more accurate method like Newton-Raphson with the added 
advantage of quadratic convergence (true quadra tic convergence 
characteristics is lost in the Decoupled method ). 

4.4 FAST DECOUPLED METHOD; 

Tables 4.5 and 4.6 give the computation time and memory 
requirement for the three systems. In this method, the 
iteration time is reduced to a great extent as compared with 
other methods. Memory requirements for this method are the 
same as that for Decoupled method. Quadratic convergence 
feature is lost in this method and thus wo need a few addi- 
tional iterations for convergence as compared to Newton-Raphsoi; 
However the increase in overall time is not much because of 
the lower per iteration time. 

4.5 Q LIMIT ADJUSTMENT: 

Q limit adjustment is also tried out using the 
following three schemes. 
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Table 4.5 


14 bus system 


; Type ; C.P.U. 

I . Time 

i i 

Iterations 

Time per 
Iteration 

Specified 

6 s 

■Achieved ■ 

■ 6 A : 

i I 

• 1 

Without ' 1.38 

Sparsity : 

|(Case I) 1 

10 

0*138 

0.001 

0.00080 1 

; 

1 » 

With i , ^ 

Sparsity , 
i(Case II) i 

j ! 

J 

10 i 

i 

1 

! 

0.122 ^ 

0.001 

j 

0.00080 1 

1 ■ 

1 j 


57 bus system 


Type 

C.P.U. Iteration® ,’Time per 

1 Specif ied 1 

i 

: Achieved! 


Time ' :Iteration 

1 6 s 

I . . ■ 

1 

iWithout i 
Sparsity j 
;(Case l) j 

9.76 i 6 : 1.626 

, i 

: ! 

1 i 

0.02 1 

1 

1 

; i 

0.019 1 

■ \ 

;V\/ith ' 

Sparsity 1 

i(Case II) i 

1 ( 

i ^ ’ 

5.28 : 6 : 0.88 1 

: : i 

CM 

O 

« 

O 

1 

{ 

0.019 


100 bus system 


■ Type 

i 

! 

1 

1 C.P.U. Iterations 
! Time 

1 

Time per 
Iteration 

: Specified 

; /■■■ 

’ - s 

Achieved; 

A ; 

1 

t 

i 

jWithout 
Sparsity 
|(Case Ij 

1 

t 

1 52.57' 

1 

1 ; 
t 

7 

7.51 

; 0.001 

1 

; i 

; i 

0.000281 

1 

1 

jvi/ith 

jSparsit'N 
(Oase II, 

i 

; 14 

1 

' i 

s 

7 

i 

2.0 

j 

; 0.001 

, 1 

0. ■00028' 

\ 

1 j 
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for Decoupled method). 
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1. P-V to P—Q switching 

2. Voltage perturbation 

3. Voltage perturbation using feedback 

The flow chart for the above mentioned methods are 
attached alongwith. 

All these have not worked out very neatly. P-V to 

P-Q switching scheme works well for the 57 bus system, 

probably because of the number of P-V buses is not very large. 

When applied to the 100 bus system with 19 P-V buses 

(including slack) violations keep occuring at every iteration 

does 

and the solution £ not converge. The addition of soft 
constraints reduced the number of violations but without 
appreciable overall gain. In the voltage perturbation method 
the voltage of the P-V bus is perturbed slightly to 0.1/^ for 
57 bus and O.b^C for 100 bus system to adjust the Q limits. 

The bus is treated as a P-V bus throughout. In this case 
convergence is obtained in an iteration when Q is being 
violated at one of the buses (bus No. 73 for the 100 bus case). 
Also, the Q at other P-V buses goes too far inside the Q 
limits. 

In the third scheme, although the Q at most of the buses 
is within the tolerance band, yet, at one of the buses it is 
completely out of limit. This is because, the solution 
converges when the Q violation takes place at one of the 
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bus. If we introduce the constraint that both should be 
satisfied simultaneously then it does not converge at all. 

All the schemes (1,2, and 3) have been tried out only in 
the case of Newton-Raphson method. 

It is clear that the success in all the' above schemes 
especially schemes 2 and 3 is due to various empirical 
adjustments. At the same time, the adjustments are system 
dependent i.e, they may work for a particular system only 
(this is true for voltage perturbation scheme no. 2). 

The flow diagrams data and load flov; results for the 
three systems by various methods have been attached 
alongwith. It is to be noted that for the 100 bus problem, 
results are for the adjusted solution with scheme (3), An 
unadjusted solution takes 4 iterations (Q’s being violated 
at 6 buses) and 16,6 seconds. The results with Decoupled and 
Fast Decoupled methods using scheme 2, are also inclosed. 

It is to be noted that the unadjusted solution will require 
lessor number of iterations. 



CHAPTER 5 


CONCLUSION 

The main objective of this thesis has been to present a 
detailed comparative study of Newton-Raphson, Decoupled and Fas1 
Decoupled methods. The importance of any load flow solution 
depends largely upon its merits regarding reliability, conver- 
gence characteristics, solution time and memory requirements. 

The above methods differ, most, in their memory and computation 
time requirements. Keeping this in view, a comparative study 
of the aforementional methods has been made. In order to 
optimize computational time, the emphasis has been on the sparsi' 
oriented programming approach. From the results obtained in 
this thesis, it is clear that this approach optimizes memory 
and/or computational time. The full potential of those methods 
is realized only when memory and computation time are optimized 
by the application of sparsity oriented programming techniques. 

For practical power system, various additional features, 
should be incorporated in the load flow program. These 
additional features are in the form of Q limits, variable 
transformer tops etc. The schemes tried out here for the Q 
limit, have not yielded satisfactory results. It is felt 
that the addition of these features in the load flow program 
would further enhance its utility. 
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APPENDIX 


The general ’continuous feedback* adjustment formula is 

X = a Ay 

where a is the ’feedback gain’ whose choice is important fc 
each type of control, each load flow method, and in some cases 
each system. The objective in choosing a is to minimise the 
total number of iterations while preserving reliable convergence. 
The slowly converging methods tend to suffer much less than the 
fast converging ones from the effects of the adjustments. The 
value of a can be choosen empirically to suit a particular 

system or else should bo approximately the sensitivity between 
X and y at the operating point. For a given system, a suitable 
fixed estimate of this can be calculated or found empirically. 

When the adjustment process is initiated, a trial 
correction (not too small) is made on the basis of an 

error more load flow iterations are then 

performed until moderate convergence is achieved, and the new 
error is y^^^. An estimate of a can now be found thus. 


a 



( sa?AB!r ) 



Fast 

Decoupled 


Flow Chart for the Hewton-Raphson, Decoupled aad Fast Decoupled 
methods using Sparsity and ordered eltoiinat aon. 
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Q-linits using error feed 'back principle wi1ii soft constraints 

NOTE’i For N bus systen. 

Let L-Q buses = MB 
and P-V Buses = HFV 
(Not including slack) 

Then after renumbering as in block two of flew chart. 
Bus No. 1 to MB will be P-Q buses 
Bus No. MB+1 to MB+NPV will be P-V buses 
Bus No. MB+liiPV+1 will be slack bus 
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